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Abstract 

o 

In a previous publication, we have constructed the Schrodinger functional in Wil- 
son's lattice QCD. It was found that the naive continuum limit leads to a well-defined 
ly-^ classical continuum theory. Starting from the latter, a formal continuum definition 

of the Schrodinger functional is given and its saddle point expansion is carried out 
to one-loop order of perturbation theory. Dimensional regularization and heat kernel 
techniques are used to determine the one-loop divergences. These are partly canceled 
Qh by the usual renormalizations of the quark mass and the coupling constant in QCD. An 

(— I additional divergence can be absorbed in a multiplicative renormalization of the quark 

boundary fields. The corresponding boundary counterterm is a local polynomial in the 
fields, so that we confirm a general expectation of Symanzik. 
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1. Introduction 



The Schrodinger functionial in a quantum field theory is defined as the transition 
amplitude between field configurations at (euclidean) times a;o = and xq = T. It can 
be written as a path integral, where only those euclidean fields are integrated over which 
assume the initial and final field configurations as boundary values. In other words, the 
Schrodinger functional can be regarded as a euclidean quantum field theory, defined 
on a space-time manifold with boundaries and (inhomogeneous) Dirichlet boundary 
conditions for the quantum fields. 

Perturbative renormalizability of a quantum field theory is usually established 
using power counting in momentum space. In the case of the Schrodinger functional, 
this is not possible because translation invariance in the time direction is violated. 
Therefore, one has to worry about the universality of the Schrodinger functional even 
in perturbation theory. At this point it is useful to recall Symanzik's work on the 
Schrodinger representation in quantum field theory [1,2]. He considered the case of 
massless scalar field theory and found that the Schrodinger functional could be 
renormalized by adding two new counterterms to the action. These are proportional 
to the local composite fields (pdocj) and integrated over the hyper-planes at a;o = 
and xq = T. 

Following Symanzik, it is plausible that this result is generic for a renormalizable 
quantum field theory. We thus expect that the Schrodinger functional in a renormal- 
izable quantum field theory is again renormalizable, after inclusion of a finite number 
of additional boundary counterterms. These are local polynomials in the fields and 
derivatives thereof, integrated over the boundary. Furthermore, they must respect the 
symmetries of the theory and have canonical dimension less than or equal to 3. In the 
SU(N) Yang-Mills theory, this expectation was confirmed to one-loop order of pertur- 
bation theory [3]. It turns out that no additional counterterm is needed in this case, 
basically because there exists no local gauge invariant composite field of dimension 3 
or less, if invariance under parity is assumed. 

It is the aim of this paper to show how renormalization works out for the Schrodin- 
ger functional in QCD. In a previous publication the Schrodinger functional has been 
defined in Wilson's lattice QCD [4]. It was shown that the naive continuum limit of 
the lattice action leads to a sensible classical continuum theory. Taking the classical 
theory as starting point, we will use dimensional regularization to analyze the QCD 
Schrodinger functional to one-loop order of perturbation theory. 

The paper is organized as follows. Sect. 2 reviews some aspects of the classical 
continuum theory. In sect. 3, a formal continuum approach to the Schrodinger functional 
prepares the ground for the saddle point expansion, which is carried out in sects. 4 
and 5. We then discuss the renormalization of the Schrodinger functional and draw our 
conclusions. Finally, three appendices have been included. Appendix A summarizes 
some of the notational conventions. In appendix B, the heat kernels associated with the 
fiuctuation operators are defined and some explicit expressions are given for the case 



1 



of vanishing background field. Appendix C supplies an asymptotic expansion which is 
needed to extract the divergent part of the quark self-energy. 



2. Classical Theory 

In ref.[4], a classical continuum action was derived from the lattice Schrodinger 
functional and some aspects of the corresponding classical theory have already been 
discussed there. It is worthwhile to have a further look at the classical theory since it 
provides the starting point for the formulation of the Schrodinger functional in dimen- 
sional regularization. In particular, we construct explicit expressions for the classical 
solutions of the Dirac equations, which will later play the role of test functions. 

The classical action 

The gauge part of the classical continuum action is given by 

Sg[A] = -^ [ dxo [ d^KtT{F^,F^,}, (2.1) 
^fl'o Jo Jo 

where Ffj^^, is the field strength tensor associated to the SU(iV) gauge field Afj^, 

— d^Ajy — di/A^ -\- [A^, Ajy\, (^•2) 

and the spatial vector components of the gauge field are required to satisfy the boundary 
conditions 

Afc(a;)U„=o = Cfc(x), Afc(a;)U„=T = CUx). fc = 1,2,3. (2.3) 
It will be assumed that there are Uf degenerate quark fiavors of mass m, with action 



Sf[A,i^,^]= [ dxo [ d'i^4>{x){l/)+m}^{x) 
Jo Jo 



^d^K [^l.{x)P_^l.{x)]^^^^ - /"d^x [^l.{x)P+^l.{x)]^^^^ 

</ 



(2.4) 



The covariant derivative is = d^ + A^ and the projectors P± = |(l±7o) are used to 
project on the Dirichlet-components of the quark fields. More precisely, the boundary 
conditions in the euclidean time direction are 

p+V(a;)L„=o = z'+W, P-V(a;)L„=T = /'-W, 

(2.5) 

V'(a;)P-Lo=o = ^-(x), V(a;)P+Lo=T = ^+(x), 
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whereas, in the spatial directions, periodic boundary conditions with period L are 
assumed for all fields. 

As pointed out in ref.[4], the form of the fermionic action follows unambiguously 
from the boundary conditions (2.5), if one assumes parity invariance of the action | 
and the existence of smooth classical solutions to the equations of motion. 

Classical solutions 

The classical solutions and are determined through the boundary quark fields. 
In order to construct them explicitly, one needs the classical quark propagator, defined 
through 

{lp+ m)S{x, x') = 8{x - x'), (2.6) 

and the boundary conditions 

P+S{x,x'%,=o = 0, P-S{x,x%,=T = 0, 

(2-7) 

S{x,x')P_l,^^ = 0, S{x,x')P+l,^^= 0. 

In ref.[4], it has been shown that this propagator is well defined if the gauge potential 
vanishes. It will be assumed here that this remains true in presence of a sufficiently 
well-behaved gauge potential. For the classical solutions one then obtains 



Mx) = J d3x'[5(a;;0,x')/>+(x') + 5(a;;T,x')/>'_(x')], 
Mx) = J d3x'[^_(x')5(0,x';a;) + ^V(x')5(T,x';a;)]. 



(2.8) 



After partial integration in eq.(2.4), the action of the classical fields assumes the more 
symmetric form, 

d'^[AiAi],^^o-2 ^'^["PciM^^^T^ (2.9) 
and, using eq.(2.8), a functional of the boundary quark fields is obtained. 



Sf[A,^,uM = - j d'xdV 



^_(x)5(0,x;T,x')/>'_(x') 

^_(x)5(0,x;0,x')/>+(x') 
^V(x)5(T,x;T,x')/>'_(x') 

^V(x)5(T,x;0,x')/>+(x') 



(2.10) 



I Parity invariance of the QCD action is assumed throughout the paper. In particular the vacuum 
angle 6 [5,6] is taken to vanish. 
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3. Formal continuum approach 



In the continuum, the Schrodinger functional can be represented as a formal eu- 
clidean path integral [3], 

Z[-p\,p'_,C'rP-,P+,C] = j D[A] j DiV'PiV'PiAle-^C^'^'^]. (3.1) 

S = Sg -\- Sf IS the euclidean action in the continuum (2.1), (2. 4) and the functional 
integral extends over those fields which satisfy the boundary conditions 

M^)\.,=, = Cfc^(x), M^)L=T = cUx), 

P+V(a;)L„=o = A(x)/>+(x), P_^x)l^^^ = />'_(x), (3.2) 

H^)P-L=o = /^-(x)A(x)-S ^l^{x)P+\^^^^ = ^V(x) . 
Here, denotes the gauge transform of the boundary gauge field C, 

= ACfcA-i + A9fcA-\ (3.3) 

with the time independent gauge function A(x). The appearance of the gauge trans- 
formed boundary fields in eq.(3.2), together with the integration over the gauge function 
A(x) corresponds to a projection on the physical subspace of gauge invariant wave func- 
tions. On the lattice, this projection is equivalent to the integration over the timelike 
link variables [4]. 

All fields are periodic in space. To preserve periodicity under a gauge transforma- 
tion, only periodic gauge functions are admitted. In particular A(x) has to be periodic 
and can thus be interpreted as a mapping from the 3 dimensional torus to the group 
SU(iV). Such functions fall into disconnected topological classes which are labeled by 
an integer winding number 

n = ^ ^ d^x eki, tr { (ASfcA"!) (AS^A-^) (AS.A-^) }. (3.4) 

It is possible to convert the integral over A to a sum over winding numbers by making 
use of gauge invariance [3]. To this end fixed representatives A„ are chosen for gauge 
functions with winding number n. The Schrodinger functional then reads, up to a field 
independent normalization factor, 

oo . 

Z[p'^,p'_,C';p.,p+,C]= /D[V']D[V']D[A]e-^[^'^'^]. (3.5) 

n= — oo 

For given n, the boundary conditions (3.2) are to be taken with the representative A„ 
instead of A and for later convenience we make the choice Aq = 1. 
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Minima of the action and gauge group 

At small values of the gauge coupling go, the dominant contributions to the functional 
integral are expected from small neighborhoods around the absolute minima of the 
action. We make the assumption that the minima of the action are determined by the 
pure gauge theory alone. In this respect quark fields only play a secondary role. 

The absolute minima of the pure gauge action depend on the choice of boundary 
gauge fields. In the following, we consider only those boundary fields for which the 
minimizing configuration is unique (up to gauge transformations) and attained in the 
winding number n = sector. The existence of such boundary fields has been estab- 
lished in ref. [3]. It is the corresponding minimal action configuration B^i^x) which is 
referred to as the induced background gauge field. 

With these assumptions, the gauge group is essentially determined as the group 
Q of gauge transformations which leave the boundary gauge fields intact. In order to 
further specify Q we assume irreducible boundary gauge fields, i.e. any gauge function 
A(x) for which = C must be constant and proportional to the identity. A SU(iV) 
gauge transformation G ^ hence satisfies 

\ Zm at Xq = 0, 
''(^^=lz„, at.o=T, ^'-'^ 

with Zm = exp2xim/iV, and some integer numbers m and m'. Since constant gauge 
transformations act trivially on gauge fields, the group which is relevant for the gauge 
fixing procedure is Q = Q /Z^. Here, denotes the center of SU(iV) which consists of 
all elements Zm, m = 0, .., N — 1. Q can thus be identified with the m' = component 
of Q. This group acts freely on gauge fields and has N disconnected components labeled 
by m. If we choose fixed representatives 0^ for each component Qm of Q, each element 
g G Qm can be uniquely represented as 

g = Qrnh, (3.7) 

where h is an element of Qo, the component containing the identity. 

Degeneracy of the saddle point 

In perturbation theory, one chooses a point on the gauge orbit of the minimal action 
configuration B and parametrizes the infinitesimal fiuctuations around this minimum, 
A = B + goq. The gauge fixing procedure then separates the infinitesimal gauge direc- 
tions from the physical field fiuctuations and hence defines a saddle point. 

Since the gauge group Q consists of N disconnected components, any gauge orbit 
decomposes accordingly. More precisely, the gauge orbit of a gauge potential A is the 
union of the orbits of all A^"^ under the identity component Qq. Hence, the gauge fixing 
procedure selects a saddle point on each component of the gauge orbit. This does not 
cause any problem in the pure gauge theory, because each component makes the same 
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contribution to the functional integral, i.e. the saddle point has an iV-fold degener- 
acy. However, when quark fields with non- vanishing boundary fields are included, the 
degeneracy is lifted, because the boundary conditions for the quark fields are not the 
same on the diff"erent components of a gauge orbit. 

A nice method to deal with this situation starts with the observation, that a gauge 
transformation with 0^ leads to 

Z[-p\,p'_,C';-p-,P+,C] = Z[p'^,p'_,C';p.z-\zmP+,C]. (3.8) 

One may take the average over all values of m, 

^ = nT. Z[p'+,p'-,C';p-z-\zmP+,C]. (3.9) 

m=0 

Next, we shift the fermionic integration variables in order to integrate over fields sat- 
isfying homogeneous boundary conditions. The dependence on the quark boundary 
fields and the boundary values Zm of the gauge functions now resides in the action and 
is displayed by a superscript {m), 

Z[p'^,p'_, C';-p-,P+,C] = j^Y. mmme-^ [A,x + v,x + v]_ (3.10) 

m=0 

If the boundary values Zm are interpreted as additional dynamical variables of the 
action, transforming as 

under a gauge transformation 0;, the enlarged action is gauge invariant, i.e. 

5(™+')[A"SxOr\0,x] = 5(™)[A,x,x]. (3.12) 

Furthermore, the summation over m can be regarded as part of the measure, so that 
the measure and the action are separately invariant under the action of the gauge group 
Q. This is the standard situation to which the gauge fixing procedure can be applied as 
usual. Leaving the sum over m until the very end, one can treat the functional integral 
for fixed value of m in exactly the same way as the pure gauge theory functional. In 
particular, one recovers the iV-fold degeneracy of the saddle point, because the fermionic 
boundary conditions now transform covariantly from one component of a gauge orbit 
to any other one. 
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Dimensional regularization 

The formulation of the dimensionally regularized theory starts with the extension of 
space-time to a Z)- dimensional manifold. For the p additional dimensions we choose 
the torus = 5^ X 5^ X • • • X 5^. Periodic boundary conditions are assumed for all 
fields so that the additional dimensions are analogous to the spatial ones. 

As usual in dimensional regularization, external momenta and fields have physical 
components only. Concerning the Schrodinger functional, this means that the boundary 
fields are independent of the extra dimensions. Consistency then requires that also the 
gauge functions G ^ and A„ depend on the physical dimensions only. 

Using the conventions of appendix A, the action in Z) = 4 + p dimensions reads 



x + v,x + v] = S,[A] + Sf[A, v, v] + ^"^[A, x, x], 



(3.13) 



with 





(3.14) 



Sr[A,x,x] 




xo = 




The fields % and % are solutions of the equations of motion which follow from S'^™^ 
One finds 




and a similar expression for %. The Z)- dimensional quark propagator 



S' 



[x, x') satisfies 




(3.16) 



For the action, one derives the expression [cf.eq.(2.10)] 



d^-^i d^-^i' ^_(x)5'(0, i; 0, i')/'+(x') 



+ ^V(x)5'(T,i;T,i')/>'_(x') 
+ ^_(x)5'(0,i;T,i')/>'_(x') 



(3.17) 



+ z„^V(x)5'(T,i;0,i')/>+(x') , 



7 



where the dependence on Zm is explicit. Taking into account the behavior of the 
propagator under a gauge transformation 0, 

S'{x, x') ^ n{x)S\x, x')n{x')-^, (3.18) 

one easily verifies that indeed both parts of the quark action in eq.(3.14) are invariant 
under the action of the gauge group Q. 

Gauge fixing 

The gauge fixing procedure may entirely be taken over from refs.[3,7]. Due to the 
assumptions on the minimal action configuration, only the winding number n = 
sector contributes to the saddle point expansion. Restricting our attention to this 
sector, the gauge field A is decomposed as follows, 

A^{x) = B^{x) + goq^{x), Ap,{x) = goqp,{x). (3.19) 

The spatial components of the fiuctuation field q hence satisfy homogeneous boundary 
conditions, 

g^(0,i) = 0, g^(T,i) = 0, (3.20) 

whereas the time component go remains unrestricted at this point. 

The so-called background gauge is implemented by choosing the gauge fixing term 
(with bare gauge fixing parameter Aq) and the Faddeev-Popov action as in ref.[3]. 
Denoting the covariant derivative in the adjoint representation by V p, = dp, + Ad Bp, 
we have 

Sgi[B,q] = -Ao / d^xtT{VpqpVi,qi,}, 

(3.21) 

SFp[B,q,c,c] = 2 / d^xti{cVp{Vp + goAd qp)c}. 
The total gauge fixed action reads 

'^tr^ai[^'?'C,c,x + t;,x + t;] =Sg[B + goq] + Sgi[B,q] + SFp[B,q,c,c] 

(3.22) 

+ Sf[B + goq,v,v] + s''P\b + (?og,x,x]. 

A careful analysis shows that the Faddeev Popov ghosts and the time component of 
the quantum field q satisfy the boundary conditions [3], 

Vogo(ic) = ci^x) = ci^x) = at a;o = and xq = T. (3.23) 
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Finally, the gauge fixed, dimensionally regularized Schrodinger functional is given by 



JV— 1 

^=0-^ (3.24) 

and provides the starting point for the saddle point expansion. 

The role of 

Before turning to the saddle point expansion let us understand the role of the variables 
Zm- The summation over m may actually be carried out if one expands the exponential 
of the z^-dependent part of the action, viz, 

exp{-4™)[5 + ffog,X,x]} = E^{4™^[^ + ^o«'^'^]}'- (3-25) 

fc=0 

In view of eq.(3.17), this corresponds to an expansion of the Schrodinger functional in 
powers of the boundary quark fields. 

The = contribution is independent of z^ so that the sum over m is trivial. For 
k = 1, one sees that due to 

JV-l 



m=0 



only the part which contains either the fields at a;o = or the fields at a;o = T survives 
the summation over m. Concerning a general term in the expansion, one still obtains 
the contributions which involve only boundary fields at one of the boundaries. In 
addition, one may have pairings z^z^^ = 1, and powers thereof. If we use the relation 

^ri, if l = nN, n=0,l,2... 
^ l^^mi \ Q otherwise, ^ ' 

m=0 ' ' 

we also see that terms which contain z^ = 1 or powers thereof contribute to the 
functional integral. 

The physical interpretation is obvious. Only gauge invariant ("colorless") states are 
propagated in euclidean time. These are states, formed from the boundary quark fields, 
which have the quantum numbers of either "mesons" or "baryons" or are combinations 
of both. This physical picture is thus seen to be a consequence of gauge invariance. 



9 



4. The saddle point expansion 



Having understood the role of the variables z^, one may leave the summation over 
m until the very end and restrict attention to the gauge fixed functional integral with 
fixed value of m, 

Z("^^[p'^,p'_,C';p-,P+,C]= [ D[g] / D[c]D[c] / DMD[t;] 



X e 



From the technical point of view, the evaluation of the Z^™) is the same for all values 
of m. In what follows we will therefore restrict attention to Z^'^\ Eventually, all other 
contributions may be obtained by appropriately replacing the boundary quark fields. 
It proves convenient to define an auxiliary quantity T through 



T[BJ,l,^l^,l] = -lnZ^'^[p'^,p'_,C';p.,p+,C]. (4.2) 

By slight abuse of language, r[5,^cbV'cz] will be called the eff"ective action in the 
following. Its arguments are the classical fields of sect. 2, with Afj^ replaced by the 
background field Bfj^. Due to the relation 

X{x) = il^ciix) - go J d^xiS'{x,xi)i{xi)il)ci{xi), (4.3) 

and an analogous one for %, the classical fields are equal to % and % for vanishing gauge 
coupling go and Zm = 1- Using eq.(4.3), one may rewrite the action as a function of 
the classical fields, viz 

sf\B + goq,X,x\ =sf\B + goq,tpci,tpci] 

+ 9o J d^x-4)ci{x)i{x)il)ci{x) ^4 4^ 

- gl d^Xld^X2^l)cl{xl)i{xl)S'{xl,X2)i{x2)^l)cl{x2). 



The interesting property of the eff"ective action is its field dependence. Henceforth, 
we will drop field independent additive contributions without further notice. This 
corresponds to the omission of field independent multiplicative contributions to Z^'^\ 
At the end, one therefore has to make sure that all components Z^™) are normalized 
in the same way. This can be achieved, e.g. by referring to their common value in the 
case of vanishing quark boundary fields. 
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Background field gauge symmetry 

The effective action is invariant under arbitrary gauge transformations of the back- 
ground field and the classical quark fields, 

r[5",V;,iO-\OVd] = T[B,4>,i,^,i]. (4.5) 

To see this explicitly, first note that the total gauge fixed action (3.22) is invariant under 
the combination of a gauge transformation for the background fields and covariant 
rotation of the fiuctuation fields, 

qij,{x) Q{x)qp,{x)Q{x)~'^, 
c{x) 0(a;)c(;c)0(a;)-\ 

c{x) 0(a;) c{x) 0(a;)-\ (4.6) 
v(x) v(x) 0""^, 
v[x) fl[x) v[x). 

The covariant rotation corresponds to a unitary transformation in the spaces of fields 
integrated over. The functional measure is left invariant and the gauge invariance of the 
eff"ective action hence follows. In particular, it should be emphasized that the eff"ective 
action is independent of the gauge fixing parameter Aq. One is therefore free to set 
Aq = 1 which is the most convenient choice for the perturbative calculations to be 
carried out later. 

Expansion in powers of go 

The eff"ective action has an expansion in powers of the coupling constant go, 

r = g^'To + Ti + glT2 + 0{gt). (4.7) 

To determine the coefficients of this expansion, we first expand the action in powers of 
the bare coupling. To 0(1), the action reads 

s['2JB, q,c,c,x+v,X + v] = Sg[B] + S^^'^B, V'ch V'd] + [ d"" x v {ip + m) v 

. (4.8) 
-2 / dPxii{\q^{k^ g)^ + cAo c] + 0{go). 

Denoting the field strength tensor associated to the background gauge field B^ by G^^^, 
the fiuctuation operators for the gluon and ghost fields are defined through 

(Aig)/i = -V^V^g/i + (1 - Ao)V/iV^g^ - 2[G/i^,g^], 

(4.9) 

Aqc = -Vi,Vi,c. 
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The first order term of the effective action is the gauge part of the classical action, 

To [B, V'd] = To [B] = gl Sg[B] (4.10) 

The next order term Ti receives contributions from the determinants of the fluctuation 
operators for ghost, gluon and quark fields. To have a unified notation, we define the 
operator A2, 

A2 ='(^+m)(-^+m), (4.11) 

with the covariant derivative Dj^ = djj,-\-Bj^. A2 is defined on spinors vi^x) which satisfy 
the boundary conditions 

P-v{^)l,=o = 0, {Do - m)P+v{x)l^^, = 0, 

(4.12) 

P+v{^)Lo=T = 0' (^0 + m)P_v{x)l^^^ = 0. 
It is not difficult to verify that the effective action to this order is formally given by 
ri[5,V;cz,Vd] = S^p[B,xl^,i,xl^,i] + ilndet Ai -Indet Aq - ^Indet A2. (4.13) 



Evaluation of the determinants 

The operators A^ [i = 0,1,2) are symmetric, elliptic (for Aq > 0) and bounded from 
below. Aq is positive and our assumptions on the background field are such that this 
is the case for Ai and A2 as well. All three operators have a complete set of smooth 
eigenfunctions and a discrete spectrum of positive eigenvalues. Their heat kernels 
(cf. appendix B) are therefore well-defined and can be used to define the logarithm of 
the determinants as meromorphic functions of p = —2e. For the pure gauge theory, 
this has been carried out in ref.[3]. Neglecting terms which vanish with e, the result is 



In det Ao 
In det Ai 



In 4x — 7e 



— h In 4x — 7e 



i-2^ao(A,)-C'(0|Ao), 



i-2^ao(Ai) + 2i-2^ao(Ao) - C'(0|Ai) 



(4.14) 



The derivative of the zeta function C(s| A) at zero is given by 

C'(0|A)=ao(A)7E+ / - Tre-*^- J]i-^/2a,/2(A) 
Jo ^ L 



^(i-l)ao(A) 



(4.15) 
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The coefficients 0^/2 (A) appear in the Seeley-DeWitt expansion for the trace of the 
heat kernels [8,9], viz 



Tr e-*^ a2(A) t-' + a3/2(A) i'^/^ + ai(A) + . . . . (4.16) 

Since A2 has the same properties as Aq, we need not discuss it separately and infer its 
contribution from eq.(4.14). In order to obtain the singular part of the effective action 
it thus remains to determine the Seeley coefficient ao for the 3 fluctuation operators. 

The Seeley coefficients can be calculated using known techniques [10,11]. They 
get two kinds of contributions. The so-called volume terms contribute only for even j 
and are given as local polynomials in the gauge fleld Bfj^ and its derivatives, integrated 
over the volume. The boundary terms may contribute to all terms except the leading 
one (j = 4). In contrast to the volume terms, the local composite flelds have support 
only at a;o = or at a;o = T and are integrated over the corresponding hyper-planes. 
On dimensional grounds, the canonical dimension of the local composite flelds must be 
4 — J for the volume terms, and 3 — j for the boundary terms. 

It is at this point that one can take full advantage of the fact that the background 
gauge does not violate the gauge symmetry in the background fleld. Therefore, the 
coefficients which appear in the Seeley-DeWitt expansion are gauge invariant. On the 
other hand, a boundary contribution to ao ( A) must arise from a dimension 3 polynomial 
in the background fleld and its derivatives. Since such a gauge invariant polynomial 
(with even parity) does not exist, one may conclude that ao only receives the usual 
contribution from the volume term. 

Using the known result (see e.g. [12]), one gets 

N f 

ao(Ao) = ^ j A^xtv{G^,G^,}, 

ao(Ai) = -20ao(Ao), (4.17) 

472. f 

ao(A2) = --^ao(Ao). 



Renormalization of the coupling constant 

Using the result (4.17), we notice that the divergent part of the eff"ective action is 
proportional to the classical action of the background gauge fleld, 

ri[5,v;,;,Vcr=° ^^4T^ro[5] + o(i). (4.18) 

The singularity will be canceled by the usual coupling constant renormalization. To 
see this, let /x be the normalization mass and deflne the renormalized coupling g-^a, in 
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the minimal scheme (MS) of dimensional regularization as usual through 



9 ms 



(4.19) 



1=1 



According to the convention of the minimal scheme, the singular coefficients contain 
only the poles in e without constant part, i.e.. 



(4.20) 



k=l 



The effective action near four dimensions as a function of the renormalized coupling is 
then given by 



r[5, V'cb V'cz] = 



-2e 



7^ - ^i(s) 



To[B] + T^[B,i,,i,i,,i] + 0{gl,s)- (4-21) 



The effective action is hence finite up to terms of the order 0{gl/ig), provided one sets 

1 llN-2nf 



^i(s) 



3£ 16x 



2 ' 



(4.22) 



which is the same coefficient as usual [13-15]. In 4 dimensions, the field dependent 
parts of the effective action to this order are then given by 



1 lliV -2nf , T iV , ^ 



D=i I fl'MS 48x2 



48x2 



1 z'' 



C'(0|Ai) + C'(0|Ao) + K'(0|A2) 



(4.23) 



5. One-loop boundary effects 



In order to observe the effect of the boundary one has to look for one-loop graphs 
involving the quark boundary fields, or, equivalently, the classical quark fields. It turns 
out that only 1 graph, the quark self-energy needs to be calculated explicitly. Its 
evaluation is done using a method which is originally due to Liischer [12] and may be 
adapted to the case at hand. 
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Preliminaries 

The effective action is a functional of the classical quark fields which play the role 
of external sources. Expanding the eff"ective action in powers of the classical fields, 
one may define n-point functions as the coefficients of this series. In order to identify 
possible boundary counterterms recall that these must be given by local composite 
fields of canonical dimension 3. Therefore, it is sufficient to consider only the 2-point 
function. The corresponding contribution to the effective action should refiect the form 
of the boundary counterterms, which are expected to cancel the boundary divergences 
in the higher n-point functions as well. 

To identify the graphs which contribute to the 2-point function, we inspect the 
interaction part of the gauge fixed action (3.22). Taking into account the couplings 
to the classical quark fields (4.4), one can construct 2 types of one-loop graphs which 
contribute to the 2-point function. First, a tadpole graph with a ghost, gluon or quark 
field in the loop is obtained, due to the coupling of the classical quark fields to a single 
gluon. We will argue that no divergent part is to be expected from the tadpole graphs. 
In order to present the argument, it is however useful to first acquire some experience. 
This will be done by considering the second contribution at this order, which is the 
self-energy graph for the quark field, 

T2[B,-4)ci,il'ci] = -gl J d^xid^X2-4)ci{xi)'S{xi,X2)il)ci{x2) + ■ ■ ■■ (5.1) 

Explicitly, the one-loop expression for the self-energy S is given by 

t{x, x') = D{x, xX^^i.T'^Sii, x')^i,T\ (5.2) 

Here, I)i^x,x')'i^^ is the Z)- dimensional free gluon propagator. It is related to the heat 
kernel through 

/>oo 

D{x,x')% = / dtkt{x,x'\kr)%. (5.3) 

To have a similar representation of the quark propagator it is convenient to define the 
integral kernel 

fCt{x, x'\k2) = m)Kt{x, x'\A2), (5.4) 

which is defined such that 

/> oo 

S{x,x')= / dtlCt{x,x'\A2). (5.5) 
Jo 

In the following, the divergent part of the quark self-energy will be evaluated at one-loop 
order. 
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The method 



The standard techniques of dimensional regularization [16-18] work in momentum space 
and thus rely on translation invariance. The latter is broken by the presence of the 
boundary as well as by the background fields. We therefore follow the method proposed 
in ref.[12], which may be adapted to the case of manifolds with a boundary. 

The procedure is briefiy described as follows. One first writes down the expression 
for a given Feynman graph in position space according to the appropriate Feynman 
rules. Next, one inserts the heat kernel representation of all propagators and performs 
a Fourier transformation in the p extra directions. Taking into account translation in- 
variance and the vanishing of external momenta in these directions, the p-dimensional 
integrals are either trivial due to ^-functions, or gaussian and can be calculated explic- 
itly. The outcome is an integral over the proper times with an integrand consisting of 
the purely 4-dimensional heat kernels, multiplied by an explicitly known, p-dependent 
function of the proper times. The latter regulates the proper time integrations at the 
lower end, for suitable choice of p = —2e, (usually for sufficiently large Ree). Smear- 
ing the resulting expression with a test function, an ordinary function of p = —2e is 
obtained. To analytically continue this function to e = 0, one must find out for which 
values of the proper times the integrand is not integrable without the regulating factor. 
This is usually the case when some or all of the proper time parameters simultaneously 
approach zero. By subtracting and adding the asymptotic expansion of the integrand 
in these regions, one may isolate the poles in e and eventually obtain the Feynman 
diagram as a meromorphic function of e. 

Application to the quark self-energy 

Using the notational conventions of appendix A, we define the Fourier transform in the 
extra dimensions. 



The self-energy is a distribution which is to be smeared with a smooth and spatially 
periodic test function f[x,x'). At the boundary, / may take any finite value and it 
is assumed that all its derivatives at the boundary exist. The test function has color, 
fiavor and spinor indices which match those of the self energy so that S(/), defined by 



is an ordinary function of p = —2e. Note that the contribution to the eff"ective ac- 
tion (5.1) has exactly this structure, provided the classical quark fields are considered 
as test functions. This will be used in sect. 6, when the renormalization of the eff"ective 
action will be discussed. 

When the heat kernel representation of the propagators (5. 3), (5. 5) is inserted in 




(5.6) 




(5.7) 
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eq.(5.2), one gets 



/> oo 

±{x,x') = / dhdt2Kt,{x,x'\A^)f^jf,T''ICt,{x,x'\A2)ji,TK (5.8) 

Using the factorization properties of the heat kernels (appendix B), the integration over 
the extra dimensions and smearing with a test function yields 

/>oo 

S(/)= / dhdt2[r{h+t2Wl{f,h,t2), (5.9) 
with the explicitly known function r[t), 

oo 2 2; 

r{t) = (4xi)-'/' ^""^ ^ ^'^^ (5-10) 

n= — OO 

and the integrand / = /i + /o? given explicitly by 

/i(/, h,t2)= [ d^x d^x' tr {f{x, x)Kt,{x, x'\A^fX^ ^T^Kt^x, x'\A2)i.T''] , 

\ (5.11) 
/o(/, h,t2)= I d^x d^x' tr {f{x, x)Kt,{x, x'\Aor%T''ICt,{x, a;'|A2)7^T^} • 

The contraction of the p additional 7-matrices in 7o(/?^i?^2) yields a factor of p, so 
that the dependence of the integrand on p = —2e is completely explicit. 

Properties of the function /(/, ti,t2) 

I[f,ti,t2) is an ordinary function of the proper times ti and t2. Due to the positivity 
of the fluctuation operators, I[f,ti,t2) vanishes exponentially for large values of one 
of its arguments. Concerning its behavior at small values of the proper times we flrst 
recall that the heat kernels satisfy 

lim J d^xd^x'ti if {x,x')Kt{x,x'\A)} = J d^xti f{x,x). (5.12) 

Note that the heat kernels themselves and their derivatives are, for non-zero proper time 
t, valid test functions. If i2 > is held flxed, we may therefore conclude that /(/, 0, t2) 
exists. Using similar techniques as in appendix C, one may also show that /(/, ii,0) 
exists for ti > 0. Hence, if one of the proper times in the function I[f,ti,t2) is kept 
difl"erent from zero, the other can be taken to zero. This yields a smooth function of 
the remaining proper time, which diverges when the latter approaches zero. Together 
with the exponential fall-ofl" at large proper times this leads us to conclude that the 
function I[f,ti,t2) is integrable everywhere, except near the origin ti = t2 = 0. 
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One-loop result 

It is convenient to perform the change of variables 



ti=tSi, t2=tS2, Si + S2 = 1, (5.13) 

so that the self-energy (5.9) reads 

S(/)= / dtt[r{t)]-^' dsids2S{si + S2 -l)I{f,tsi,ts2). (5.14) 
Jo Jo 

In appendix C, the asymptotic small t behavior of I[f,tsi,ts2) is determined, 

I{f,tsuts2) A(/,si,S2)i-' + 0(i-^/2), (5.15) 

and an explicit expression for the coefficient function A[f,si,S2) is obtained. 

Since the regulating function r(t) (5.10) behaves asymptotically as t~^^^ , the whole 
integrand in eq.(5.14) behaves as t^~^ at the lower end of the t integration. It hence 
follows that the self-energy (5.14) is a priori only defined for Ree > 0. In order to 
analytically continue to e = 0, we first split up the range of the i-integration at i = 1. 
The integral for i > 1 is analytic in e. In the remaining integral we subtract and add 
the asymptotic expansion (5.15) of the integrand. The integral over the subtracted 
integrand will exist for Ree > —1/2 and the integral over the asymptotic expansion 
can be calculated explicitly. Near 4 dimensions, we obtain 



1 

S(/)^=° (-+ln4x) / dsA{f,s,l-s) 
^ Jo 

fiOO 

+ / dhdt2[l{f,h,t2)-e{l-h-t2)A{f,h,t2)]+0{e), 
Jo 



(5.16) 



where we have used that A[f, ii , i2 ) is a homogeneous function of degree — 2 in ii and t2 ■ 
We will be mainly interested in the divergent part of the self-energy which is explicitly 
given by 



-4m J d'^xtr f[x,x) 



J d^xtvi^^^^^ f{x,x) - ^^^'^ f{x,x) - ^{x)f{x,x)} (5.17) 
I J d3xtr{/(0,x,0,x)+/(T,x,T,x)} +0(1). 
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Tadpole graphs 

Before turning to the renormalization of the effective action and the Schrodinger func- 
tional, we come back to the tadpole graphs mentioned at the beginning of this section. 
In these graphs, a gluon line connects the classical quark fields with a loop, in which 
ghost, gluon or quark fields may circulate. The contribution of all tadpole graphs to 
the eff"ective actions can be parameterized by a function Tfj^ 

r2[5,v;chV'd]L,dpoies=^"''/ ^^xi>cnMci, (s.is) 

which may depend on the background gauge field B and the quark mass m. Fur- 
thermore, Tfj^ must be gauge covariant, due to the background gauge symmetry of the 
eff"ective action (cf. sect. 4). Note also that Tfj^ need not transform as a Lorentz vector, 
because Lorentz invariance is violated by the presence of the boundary. Fortunately, 
it is not necessary to treat the tadpole graphs explicitly. For, we are interested in the 
divergent contributions to the eff"ective action. Such contributions cannot arise from 
the tadpole graphs as the following argument shows. 

The lesson that can be learned from the treatment of the self-energy graph is, that 
only local terms appear in the asymptotic expansion. Possible divergences may there- 
fore be identified using naive power counting. When applied to the tadpole graphs, 
this means that a divergent contribution to the function Tfj^ must have canonical di- 
mension 1. The requirement of gauge covariance for Tfj^ excludes terms which involve 
the background gauge field. There could be a term proportional to the mass or a 
^-function, confining the term (5.18) at the boundary. 

At this point, the crucial observation is that, for vanishing background field B, 
the tadpole graphs vanish identically. The reason is that each vertex contains a SU(iV) 
generator and all propagators are, for vanishing background field, diagonal in color 
space. The tadpole graphs therefore become proportional to the trace of a SU(iV) 
generator in either the fundamental or the adjoint representation (see appendix A). 
Consequently, the tadpole graphs can only contribute finite terms to the eff"ective action 
and may therefore be ignored as far as renormalization is concerned. 
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6. Renormalization 



Renormalized effective action 

All divergent one-loop contributions to the 2-point function are assembled in the self- 
energy (5.16). To find the corresponding contributions to the eff"ective action, we note 
that the classical solutions ipd and ipd have all the properties which are required for a 
test function (cf. sect. 5). We define the special test function, 

fcl{x, x')st;ij;al3 '= i'cl{x)sia'>Pcl{x')tjl3- (6.1) 

Recalling eqs.(5.1),(5.6) and (5.7), the contribution to the eff"ective action is given by 

T2[B,4>d,^ci] = -i-'^S(/cO + ---- (6.2) 



If the test function fd is inserted in the pole part of the self-energy (5.17), one obtains 

^0 Cp _ f (j4^ ^^^^^^^^^^^^ 

(6.3) 



s(/cO' 



16x2£ 



d^x{ [■4)ci{x)il)ci{x)] 



xo = 



[^cl{x)^l'cl{x)]^^^^} 



0(1). 



To arrive at this equation, we have made use of the equations of motion for the classical 
quark fields. The pole contribution to the eff"ective action is expected to refiect the form 
of the counterterms. The first pole has the form of a mass term and is indeed canceled 
by the quark mass renormalization. To see this, we introduce a renormalized mass 
through 



n=l 



(6.4) 



and choose the renormalization constant Zm according to the prescription of the mini- 
mal scheme. Using eq.(6.4) to replace the bare mass in the action of the classical quark 
fields leads to 



fl'Ms"^j?^m^ / d^xtl}d{x)tl}ci{x) + 0{glj^^). 



(6.5) 



Since the action of the classical fields is also part of the eff"ective action [see eq.(4.23)], 
the quark mass renormalization eliminates the pole proportional to the quark mass in 
eq.(6.3), provided we choose 

^<:'(.) ^ (6.6) 
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This is the usual result for the renormalization of the quark mass in QCD [19,20]. 

The remaining divergence in eq.(6.3) is the boundary effect we have been looking 
for. One recognizes that it is proportional to the action of the classical quark fields 
[cf. eq.(2.9)]. Since the latter is bilinear in the classical quark fields, it is possible to 
eliminate the divergence by rescaling the classical quark fields appropriately. We define 
a further renormalization constant Z},, 



E^r(^)^Ms- (6.7) 

n=l 



and set 



4"(^)=^. (6.8) 



With this definition, V[B, ^^'^■ipci, Z^ ^^^V'cz] is finite to order g^^S' provided the cou- 
pling constant and the mass are renormalized in the usual way. 

Note that the renormalization constant Z}, has nothing to do with the wave func- 
tion renormalization of the quark fiuctuation fields v and v. The latter is not needed at 
all, because it would just correspond to a rescaling of the integration variables of the 
Schrodinger functional. Note however, that the usual one-loop wave function renormal- 
ization constant (for Aq = 1) is contained in the pole part of the self-energy (5.17), and 
would have played a role if we had considered correlation functions such as the quark 
propagator. 

Back to the Schrodinger functional 

We are now prepared to discuss the renormalization of the Schrodinger functional. First 
recall that the eff"ective action has been defined for the m = part of the Schrodinger 
functional (4.2). To obtain the other contributions one first expresses the classical 
quark fields through the quark boundary fields (2.8). Next, one goes through all values 
for m and multiplies the quark boundary fields with the appropriate factor z^. This 
procedure guarantees that all parts Z^™) of the Schrodinger functional are normalized 
in the same way. The physical implication of the summation over m has already been 
discussed in sect. 3. Here it is important to notice that the relation between the classical 
quark fields and the quark boundary fields is linear, for all values of m. If we introduce 
the renormalized quark boundary fields, 

^1/2 R / „l/2,R 

P+ = Z^ p^, P- = Z^ p _, 

P- = P-^ P+ = P +' 

we may therefore conclude that the Schrodinger functional, 

Z[p'^,p'^,C';-p^,P^,C] (6.10) 
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is finite to one-loop order of perturbation theory, provided the mass and the coupling 
constant are renormalized in the usual way. 



7. Conclusions 

At one-loop order, the renormalization of the QCD Schrodinger functional works 
out according to Symanzik's expectation. In addition to the usual counterterms for 
quark mass and coupling constant renormalization, a new counterterm is needed to ac- 
count for the boundary divergences. Its form coincides with the boundary terms already 
present in the classical action (cf. sect. 2). Therefore, it should not come as a surprise 
that the boundary divergence can be absorbed in a multiplicative renormalization of 
the quark boundary fields. In particular, this means that for vanishing quark boundary 
fields, the Schrodinger functional is finite after the usual QCD renormalizations. 

In the recent literature, the Schrodinger functional has been used for a systematic 
investigation of the continuum limit in the SU(2) and SU(3) Yang-Mills theory [21]. 
Since our result gives confidence that the Schrodinger functional is a universal ampli- 
tude, one may hope to extend these very successful studies to QCD. 

I would like to thank M. Liischer for numerous helpful discussions, and P. Weisz for 
reading the manuscript. 



Appendix A 



Conventions for dimensional regularization 

Lorentz indices are denoted by greek letters /x, z/, . . . and range from to 3. They must 
be distinguished from the Z)- dimensional indices /i, z>, . . . which run from to Z) — 1. We 
will also use fi,i/, . . ., which take the values 4, . . . , D — 1. For the spatial directions we 
use latin indices k,l, . . . which range from 1 to 3 and their Z)- dimensional counterparts 
k, l, . . .. We will make the convention that repeated indices are summed over, if not 
otherwise stated. 

The D-dimensional Euclidean space-time manifold is taken to be a cylinder [0, T] X 
ji3 jip ^j^gpg p = D — 4 is the number of extra dimensions. A point on this manifold 
is denoted by x with components Xfi which are further specified as follows 

■^jLi — Xp_ — Up^. (A.l) 
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The letter y will be reserved for points on the additional p- dimensional manifold 
throughout the paper. Occasionally, we also use the spatial vector x and its Z) — 1 
dimensional counterpart ic. 

f -matrices 

The 7 matrices are hermitian and defined by 

{TA'T^} = 7J = 7a' /x,i> = 0,...,i? - 1. (A.2) 

The anti-commutation relation (A.2) is all that is needed to carry out the Z)- dimensional 
7-algebra. For any Lorentz vector we occasionally use the notation 

In particular we make the convention 

I = ludf,, ^ = 7^9^, (A.4) 
where the definition dp, = 8/ dip, has been used. 
su(iV) conventions 

The Lie algebra su(iV) can be identified with the linear space of of all complex N X N- 
matrices Xij, i,j = 1, . . .N, which are anti-hermitean and traceless. In other words, 
the matrix X satisfies 

= -X, tr{X} = 0. (A.5) 

As a vector space over H, su(iV) has dimension iV^ — 1. An inner product is introduced 
through 

(X,F) = -2tr{XF}, (A.6) 

and may be used to define an orthonormal basis T'^, a = 1, . . . , iV^ — 1. Any element 
X G su(iV) can be decomposed over this basis 

X = X''T'', X'' = {X,T''), (A.7) 

with real coefficients X'^. Due to the linear structure of the Lie algebra, it is sufficient 
to know the Lie bracket between any two elements of the basis 

[T'',T^] = f'^'T', (A.8) 

in order to know them for any two elements of the Lie algebra. The structure constants 
f'^^'^ are real and totally anti- symmetric under permutations of the indices. They obey 
the relation 
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which is a consequence of the completeness of the basis. 

A representation of the Lie algebra is a vector space homomorphism R which 
preserves the Lie bracket (A. 8), i.e. 

[R{T''),R{T^)] = f^'RiT'). (A.IO) 

In the context of this paper only the fundamental and the adjoint representation ap- 
pear. For the fundamental representation one simply has R{X) = X while the adjoint 
representation is defined by R{X) = AdX. The representation space for the adjoint 
representation is the Lie algebra itself, and its action on an element Y G su(iV) is the 
commutator AdX(F) = [X, F]. The representation matrices of the basis elements T'^ 
are then simply given in terms of the structure constants, i.e. we have 

{AdT^y^ = -f'''. (A.ll) 

Finally, using eqs.(A.8),(A.9) and the normalization of the basis elements, one obtains 

j,aj,bj,a ^ T^T" = -Cf, (A.12) 

where Cp = (iV^ — 1)/2N is a constant which frequently appears in the loop expansion. 
Fields 

To enhance the readability, indices are often suppressed. However, a quark field '>p{x)sioi 
has fiavor, color and spinor indices. Flavor indices are taken from the middle of the 
alphabet ...) and range from 1 to Uf. Spinor indices (a,/?,...) take the values 
1, . . .,4. Concerning the color structure, a quark field is a vector in the fundamental 
representation, with color index i running from 1 to N . 

The gluon field, its field strength tensor and the ghost fields are elements of the 
Lie algebra su(iV) and can thus be decomposed according to eq.(A.7), 

= B^T", G^^ = Gl^T", c = c^T", c = c^T". (A.13) 

Finally, we define the fields in the adjoint representation, 

B{xr^ = -r'^B{x)% G{xrji = -r'^Gixy^,. (a.m) 
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Appendix B 



This appendix introduces the Z)- dimensional heat kernels for the fluctuation oper- 
ators of ghost, gluon and quark flelds. In the case of vanishing background gauge fleld, 
some explicit expressions are obtained. 

The heat kernels of the fluctuation operators 

The heat kernels of the Z)- dimensional fluctuation operators A^, i = 0,1, 2, are deflned 
as follows 

{e-''^^°c){xy = J d^x' Kt{x,x'\Aoy^c{x')\ 

{e-^^'q){x)l = J d^x'K,{x,x'\A^)f,q{x')l, (B.l) 

{e~^^^v){x)sia = J d^x' Kt{x,x'\A2)si<x;tjf3v{x')tjf3. 

They satisfy the diflusion equation 

{- + A)Kt{x,x'\A) = 0, (B.2) 

and from the general theory [10] it can be inferred that they are smooth functions of t, 
X and x', as long as the proper time t is positive. The limit t ^ yields, 

Vim Kt{x, x'\A) = S{x - x'). (B.3) 

The heat kernels have simple factorization properties which allow to separate the de- 
pendence on the p extra dimensions. With the convention 

A, = A,|^^^, i = 0,l,2, (B.4) 

one flnds, for Aq and A2, 

Kt{x,x'\Aor'' = Kr{y,y')Kt{x,x'\Aor\ 

Kt{x,x'\A2)si<x;tjf3 = Kj\y,y')Kt{x,x'\A2)sioi;tjl3, 



(B.5) 



The kernel Kf^(y,y') is deflned in eq.(B.8). Assuming Aq = 1, the factorization for Ai 
looks as follows, 

Kt{x,x'\A^rJl = Kr{y,y')Kt{x,x'\A^rJl, 

(B.6) 

Kt{x,x'\Ar)f^ = Kr{y,y')Kt{x,x'\Aor'S-^^. 
All other components vanish identically. 
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One- dimensional heat kernels 

Consider the one- dimensional laplacian A'^ on the circle with circumference L. It 
may be represented by the second derivative operator —d^/dx^, acting on i-periodic 
functions of the real variable x. Its heat kernel is given as the integral kernel of 



exp -tA^^ , 



Kt{x,x'\A'') = {4Tt)-'/' e-(^-^' + "^) (B.7) 



Using the notation of appendix A, the heat kernels for the laplacian on the tori T'^ and 
TP are then given by 

3 

iff (x,x')= n^*K'4iA^'), 

fc=l 

(B.8) 

4+p 



A2=4 



Next, consider the laplacian on the interval [0,T], acting on functions which satisfy 
Dirichlet (D) or Neumann (N) boundary conditions. We obtain the following one- 
dimensional heat kernels, 

Kf^'ix^x') = (4xi)-'/' {e-47(---'+2^")' _ e" ilC^+^'+^T")' }, 

n= — oo 
oo 

Kf''{x,x') = {4Tt)-'/' Y {e-4l(-°-^o+^^")Ve-4l(-+^'+2^")'}, 



n= — OO 
oo 



(B.9) 



n=0 



Kf''{x,x') = Kr{T-x,T-x'). 

Here, the first superscript refers to the boundary condition at xq = 0, the second to 
xq = T. We will also need the heat kernels with modified Neumann conditions, 

Kf^{x, a;')U=o = 0, (9, + m) Kf^{x, a;')U=T = 0, 

(B.IO) 

{d^-m)Kf''{x,x%=o = 0, Kf''{x,x')\^=T = 0. 

I did not succeed in obtaining explicit expressions for these heat kernels. However, if 
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T is taken to infinity, one finds, 



lim Kf^{x, x') = {4Tt)-'/' [e-(^ " ) /4i _ e-(^ + ^ ) /4i] 
lim Kf'ixo, x'o) = (4xi)-'/' [e-(^ " ^')V4i (B.ll) 
+ e-(^ + ^')V4i|i _ X + x')}] . 

The parameter integral F[t, x), 

F{t, x) =' r ds e-^i'' + 2"^^) - (B.12) 
Jo 

cannot be further simplified. For, using e.g. formula 3.322 in ref.[22], one gets 

F{t, x) = . /^e-(2"^* + ^)V4i [i _ <p(^mt±^)], (B.13) 

y 4:t 2vi 

where 



2 r -t' 



(l){x) = —= e"^ At (B.14) 
V"^ Jo 

is the error function. 
Vanishing background field 

The one- dimensional kernels of the previous paragraph provide the building blocks for 
the 4 dimensional heat kernels, in the case of vanishing background field. We define 

Kt{x,x'\l^oY''\^^^ = 8''''Ht{x,x% 

Kt{x,x'\/^^)f,\^^^ = S'''Ht{x,x')^,, (B.15) 

ICt{x,x'\A2)sia;tjl3\s = o = ^ sthjlCt{x, x')al3, 

with the r.h.s. further specified as follows 

Ht{x,x')f,^ = irf'(x,x') [Sf,oSiyoKf^{xo,XQ) + Sf,kSiykKf^{xo,XQ)], 

Ht{x,x') = Kr{^,^')Kf^{xo,x'o), 



ICt{x, x%p = (-^ + m)««, e-*"" iff (x, x') 



(B.16) 



{P.U.pKf''{xo,x',) 
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Appendix C 



The asymptotic small t expansion of the function I[f,tsi,ts2) [eqs.(5.9),(5.11)] is 
obtained by expanding the heat kernels in powers of the background field. 

Expansion in powers of the background field 

The free quark propagator in the presence of the smooth background gauge field B^i^x) 
satisfies 

S{x,x') = So{x,x') — J d'^xiS{x,xi)^{xi)So{xi,x'), (CI) 

where So{x, x') denotes the free quark propagator for vanishing background gauge field. 
When this equation is iterated, 

oo 

S{x, x') = So{x, x') + 5(")(a;, x'), (C.2) 



n=l 



the n-th order term is given by 

S^^\x,x') = {-ir J d'x,d'x2---d'x^So{x,x,) ^^^^ 
X ^(a;i)S'o(a;i, a;2)-^(a;2) • • ■^{xn)So{xn, x'). 

This expansion of the propagator is easily turned into an expansion of the corresponding 
heat kernel. With the notation 



/> oo 

S^''\x,x')= / dtlci''\x,x'\A2), 
Jo 

j* oo 

50(2;, x') = / dtK.t{x, x'), 
Jo 



(C.4) 



we obtain for the n-th order term 

n+l 

S, — 1) / d'^Xid'^X2 ■ ■ ■ d'^Xn 

(C.5) 

X ICts^{x,Xi)^{xi)ICts2{xi,X2)^{x2) ■ ■ ■${Xn)K.ts^^AXn,x'). 

Proceeding similarly for the gluon propagator, we obtain the corresponding expansion 



/cl"^(a;,a;'|A2) = (-i)" / dsi ■ ■ ■ ds^+i Sj - l) / d^a;id^a;2 • • • d^iCr 

Jo _i J 



mti,t2)= ^^™'"^(/,ii,i2), (C.6) 



m,n=0 



where the superscript (m, n) stands for m insertions in the gluon kernel and n insertions 
in the quark kernel. 
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Vanishing background field 

We start with the evaluation of the lowest order terms. Leaving the 7 algebra D- 
dimensional, contracting the color indices (appendix A), and using the notation of 
appendix B for the heat kernels, the integrand is given by 

I^'''\f,h,t2) = -Cf J d'xd'x'ti if {x,x')Ht,{x,x')j,i,'fj,ICt,{x,x')'fi,}. (C.7) 

As explained in appendix B, the kernel for the quark field is not known explicitly, 
even for vanishing background field. However, if the mass is set to zero, an explicit 
expression of the kernel can be obtained. Hence, for zero mass, the result to order 
reads 



5i(l -e) J d^x 



xtr{f^^f{x,x)-f^^f{x,x)] 

f (C.8) 
+ 5i(l-£) / d3xtr{/(0,x,0,x)+/(T,x,T,x)} 

-2(51 +52)(l-£/2) y"d3xtr{/(0,x,0,x) + /(T,x,T,x)} . 

To determine the infiuence of the mass we now assume the test function to have support 
only near the boundary at xq = 0. We may let T become infinite without changing 
the asymptotic behavior near xq = 0. In this limit, the heat kernel representation of 
the quark propagator is again known (appendix B). It turns out that there is only 1 
additional contribution at the order , so that the result for the massless case (C.8) 
is to be replaced by 

Cf m(2s - 4) / ,4 , 
I&'kH^ (si + S2) J 



Effect of the background field 

In order to identify the potentially divergent diagrams, it is useful to follow Symanzik [1] 
and split the propagators and the corresponding heat kernels in a free space (/) and 
a surface (s) part. By definition, the free space propagator is the propagator which 
would be obtained if periodic boundary conditions had been imposed in all directions. 
It is singular if its arguments approach each other, while the surface propagator has a 
singularity if its arguments are close to each other and to one of the boundaries. One 
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may set up a systematic power counting [23], and a moment of thought reveals the 
asymptotic behavior of I^'^'"'\f,tsi,ts2), 

(C.IO) 

Here, (/) denotes the contribution with free space kernels only, while (s) applies to 
diagrams with at least 1 surface kernel. 

We see here that indeed very few diagrams may give contributions at order . 
Besides the case of vanishing background fields which has already been treated, there 
are only 2 further contributions involving the background field. Furthermore these are 
the same as in free space and thus cannot contribute to the boundary divergences. 

It remains to calculate the 2 contributions for a single insertion in either the gluon 
or the quark kernel. Straightforward calculation leads to 



(C.ll) 



I^^'°\Us,,ts2) -^^^/d^.tr {/(.,. m.)}, 

where terms of the order t~'^^^ have been neglected. 
Final result 

When the equations (C.8),(C.9) and (C.ll) are combined, the asymptotic expansion 
for I[f,ti,t2) itself is obtained. The result can be written in the form 

I{f,tsuts2) A(/,si,S2)i-2 + 0(i-3/2), (C.12) 
where the coefficient function A[f,si,S2) is known explicitly and obtained through 

A{f,SuS2) = lim{t' [l^'''\f,tSutS2) + I^'''\f,tSutS2) 

(C.13) 

+ I^°''\f,tSutS2)]}. 
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